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ESTIMATION  OF  THEIR  ANGULAR  COORDINATES  BY 
ONBOARD  ANTENNA  SYSTEMS 

N.V.  Barkhudaryan,  A.Z.  Sazonov,  O.I.  Sukharevsky 

Kharkov  Military  University,  Kharkov,  Ukraine 
E-mail : sukharevsky @euro  .dinos . net 

The  technique  of  calculation  of  near-zone  electromagnetic  fields  scattered  by  complex 
shape  airborne  objects  is  present.  It  is  shown  that  the  direction  line  defined  by  a 
direction  finder  differs  from  true  one  and  depends  on  geometry  of  the  object,  its 
electrical  sizes  and  also  on  mutual  location  of  the  transmitting  antenna,  object  and 
receiving  antenna.  The  results  of  numerical  calculations  for  the  object  such  as  "airplane" 
are  represented.  Calculation  was  carried  out  in  a centimeter  band,  separately  in  a plane 
of  course  and  pitches  one,  radiation  - monochromatic.  The  angle  of  object  elevation  is 
constant  and  it  is  equal  to  three  degrees. 

Most  of  modem  airborne  objects  have  complex  geometrical  shape  of  its  surface.  The 
character  of  reflections  depends  from  orientation  of  object  with  reference  to  the 
direction  of  sounding.  The  scattered  field  incident  on  the  antenna  aperture  is  the  result 
of  interference  of  waves  reflected  from  separate  units  of  object  surface.  The 
contribution  to  the  resulting  field  of  components  with  different  amplitudes  and  phases 
results  in  the  not  plane  phase  front  of  the  scattered  electromagnetic  wave.  The  distortion 

of  phase  front  causes  to  deflection  of 
the  measured  target  angular  position 
from  true  one  as  the  direction  line  is 
defined  as  a normal  to  a phase  front  [1]. 
The  calculation  technique  includes 
simulation  of  the  object  surface  by  the 
system  of  triangular  facets  [2],  In  a 
Figure  1 the  model  of  a standard 
airplane  with  a wings  span  of  20  meters 
is  represented.  The  calculations  of 
scattered  field  consist  of  a numerical 
integration  of  surface  current  densities 
for  each  facet  in  barycentric 
coordinates.  The  calculations  were 
carried  out  by  means  of  special 
cubature  formulas  permitting  to 
evaluate  integrals  of  high-oscillatory 
functions  [2], 

The  electrical  dimensions  of  the  object 
are  large.  In  this  case  acceptable 
Kirchhoff  approximation.  Let  a plane  electromagnetic  wave  incidents  on  perfectly 
conducting  scatterer  with  surface  S located  in  the  free  space 
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exp(-jk(R°x)),  (1) 

where  R°  is  the  unit  vector  of  sounding  direction;  p is  the  vector  of  polarization;  k0  is 
the  wave  number  of  free  space;  £0,p0  are  permittivity  and  permeability  of  free  space 
accordingly;  x is  the  position  vector  of  object  surface  point  (Figure  1). 

If  E(x),  H(x)  are  vectors  of  total  field  then  scattered  field  will  be  defined  by  formula 

Es(x)=  E(x)-  E°(x).  Let  in  some  point  x0  somewhere  outside  of  S the  electric  dipole 
with  the  vector-moment  p is  located.  Vector  p is  arbitrary  on  the  value  and  direction. 
This  dipole  creates  a field.  Applying  a Lorentz  lemma,  we  have 

j u 3 (* o ) = j S' 1 (x )E  cl  (x/x o , p)- 1 r (x )H  oe  x (x/x  o , p}is , (2) 

where 

rtx-hxfl  ET=E-n(nE)  Ef  = E0-hfnE)  H cl  = n x H - (3) 


the  tangential  components  of  electrical  and  magnetic  fields,  ri  - is  the  unit  normal 
vector  to  the  object  surface.  For  the  perfectly  conducting  scatterer  E1  |s=0  and 

consequently  the  integrand  function  is  zero. 

The  field  of  an  electric  dipole  can  be  represented  in  the  form 

3, eT  (x  / x o , p) = — [v(pv)G  + k Ipo],  (4) 

eo 

jk  r / , . 

where  G = e 0 /4ra-;  r°  =^x  -xj/r;  r=|x  -x|  - the  distance  between  object  and 

receiving  antenna;  p^p^+p1  are  longitudinal  and  tangential  component  of 
polarization  vector.  Taking  into  account  that  in  Kirchhoff  approximation  formula  (2)  for 
a scattered  field 


where  Ss  is  "illuminated"  part  of  object  surface.  Using  the  representation  (5) 

calculations  of  scattered  field  were  carry  out  for  p = 0...900  and  £ = 75°, 85°, 95°  . For 

example,  in  Fig.  2 results  of  calculation  for  e = 75°  and  r = 200m  are  presented. 
Direction  finding  accuracy  depends  on  the  type  of  angle  discriminator  and  amplitude- 
phase  distribution  (APD)  in  antenna  aperture.  In  this  paper  we  simulated  combined 
calculation  accuracy  is  given  by  the  antenna:  in  p -plane  - amplitude  method  of 
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r(p)  = Re[usu;  exp(jTt/2)]/UsU; 

Ud  =U,  -U2 . For  example  dependences 
Figures  3,  4. 


direction  finding;  in  s -plane  - phase  one. 
Such  finding  direction  functioning  is 
described  by  expression 

U,,2(x0)=  ^s(x0)A12(x0)dSap, 

*ap 

where  A,  2 (x0 ) is  APD  in  antenna  aperture 
for  two  channels  of  receiving  antenna,  Sap 

is  aperture  area.  Deviation  of  measured 
object  bearing  from  true  one  is  determined 
by  expressions  r(c)  = Re[lJsUd  /USU*  - 
for  course-plane  and 

for  pitch-plane,  where  US=U,+U2, 
of  r(c)  andr(p)  from  angle  P are  shown  in 


The  analysis  of  Figures  3,  4 have  shown  that  in  the  course-plane  the  deviation  can  be 
obtained  the  values  1, 7-2,3  degrees.  In  the  pitch-plane  it  possible  to  obtain  15  degrees. 
For  azimuth  60°  in  the  s -plane  the  angular  deviation  is  more  than  object  angular  size. 
By  using  proposed  technique  one  can  efficiently  calculate  the  errors  in  airborne  object 
direction  finding  appearing  in  short-range  radar. 
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